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ABSTRACT 

Velocity  analysis  by  normal  moveout  encounters  problems  in  handling  dipping 
reflectors  or  lateral  variation  of  velocity.  Prestack  migration  provides  a  powerful 
tool  to  do  velocity  analysis,  which  is  based  on  the  following  principle:  the  imaged 
depths  at  a  common  location  are  independent  of  source-receiver  offset  when  the 
correct  velocity  is  used.  Conventional  approaches,  such  as  depth  focusing  analysis, 
generally  involve  iteration,  which  requires  repeated  prestack  migration.  In  this 
paper,  a  residual  moveout  method  for  velocity  analysis  on  multi-offset  data  is 
presented  that  needs  only  a  single  prestack  migration.  A  number  of  theoretical 
problems  in  this  method  are  studied.  When  the  velocity  has  a  lateral  anomaly, 
we  derive  a  formula  to  calculate  the  interval  velocity  from  the  stacking  velocity  by 
perturbation  theory.  A  suggested  data  processing  technique  based  on  our  method 
is  composed  of  prestack  migration  with  a  constant  velocity,  velocity  analysis, 
residual  moveout,  stacking,  velocity  conversion,  and  poststack  residual  migration. 


INTRODUCTION 

Normal  moveout,  typically  used  to  do  velocity  analysis  in  seismic  data  processing, 
is  robust  when  reflectors  are  flat  and  velocity  is  laterally  invariant.  However,  this 
method  encounters  difficulty  when  reflectors  are  dipping  or  velocity  varies  laterally. 
Some  geophysicists,  therefore,  have  concentrated  on  doing  velocity  analysis  by  migra¬ 
tion,  such  as  the  focusing  analysis  and  common  location  imaging  methods  (Jeannot, 
1986;  Al-Yahya,  1989).  When  the  background  velocity  is  correct,  the  imaged  depths 
at  a  common  location  are  independent  of  source- receiver  offset.  Otherwise,  if  an  in¬ 
correct  background  velocity  is  used,  the  imaged  depths  at  a  common  location  change 
with  offset.  In  this  situation,  a  residual  moveout  is  observed  in  common-location 
images  of  migrated  data.  The  principle  of  velocity  analysis  is  to  choose  a  background 
velocity  so  that  the  common- location  images  are  close  to  a  horizontal  alignment. 
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Conventional  methods  use  iteration  to  correct  velocities.  In  each  step  of  iteration, 
prestack  migration  is  required.  Consequently,  the  total  computation  is  so  huge  that 
.  application  of  this  method  is  limited. 

It  would  be  desirable  to  estimate  velocity  directly  from  the  residual  moveout.  If 
this  is  successful,  only  one  prestack  migration  is  required.  To  make  this  estimate,  we 
need  a  quantitative  relationship  between  the  residual  moveout  and  the  error  in  back¬ 
ground  velocity.  Some  geophysicists  applied  this  idea  to  velocity  analysis  (Doherty, 
1976;  Deregowski,  1990).  However,  no  general  formula  was  derived. 

In  this  paper,  the  residual  moveout  is  used  to  do  velocity  analysis  on  multi-offset 
migrated  data.  By  means  of  the  reflector  equation,  we  derived  a  general  formula  for 
the  residual  moveout.  Theoretically,  an  arbitrary  background  velocity  can  be  used  to 
do  velocity  analysis.  However,  we  prefer  using  a  constant  background  velocity  so  that 
the  residual  moveout  is  an  explicit  function  of  the  velocity  error.  Furthermore,  we  may 
use  Stolt  migration  for  constant  velocity;  this  is  much  faster  than  other  algorithms. 
Under  the  assumption  of  a  small  offset,  we  obtain  an  analytic  expression  for  the 
residual  moveout  that  is  an  explicit  function  of  background  velocity,  curvature  of 
unmigrated  data,  and  slope  of  unmigrated  data.  Furthermore,  this  residual  moveout 
is  independent  of  the  dip  of  the  reflectors  when  the  true  velocity  is  a  constant.  Parallel 
to  NMO  velocity  analysis,  we  can  use  a  semblance  and  velocity  scans  to  do  RMO 
velocity  analysis.  Therefore,  the  program  for  RMO  velocity  analysis  is  similar  to  the 
existing  program  for  NMO  velocity  analysis,  except  in  the  residual  term. 

After  velocity  analysis,  we  obtain  the  stacking  velocity  that  is  assumed  to  equal 
RMS  velocity  for  a  laterally  invariant  medium.  However,  when  velocity  has  an 
anomaly,  this  assumption  is  invalid.  Lynn  and  Claebout(1982)  proposed  a  formula 
that  gives  a  relationship  between  the  anomaly  of  the  stacking  velocity  and  the  anomaly 
of  the  interval  velocity.  But  this  formula  requires  a  stratified  medium.  In  this  pa¬ 
per,  we  derive  a  formula  that  is  valid  for  arbitrary  velocity.  This  formula  computes 
the  anomaly  of  the  interval  velocity  from  the  anomaly  of  the  stacking  velocity  by  an 
integral  equation. 

We  have  two  ways  to  do  migration:  residual  prestack  migration  or  residual  post 
stack  migration.  We  prefer  the  latter  because  poststack  migration  runs  much  faster 
than  prestack  migration.  The  residual  moveout,  stacking,  and  residual  poststack 
migration  are  recommended  to  obtain  the  final  structural  image. 

REFLECTION  EQUATION 

Seismic  signals  consist  of  amplitude  information  and  phase  information.  Under  the 
high  frequency  assumption,  phase  information  is  simplified  to  traveltime  information. 
Therefore,  seismic  signals  can  be  described  approximately  by  traveltime-offset  curves. 
A  traveltime-offset  curve  in  seismic  record  corresponds  to  a  specific  reflector  in  the 
earth.  Given  a  reflector,  we  can  compute  the  traveltime-offset  curve  by  modeling. 
On  the  other  hand,  given  a  traveltime-offset  curve,  we  also  can  compute  the  reflector 
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position  (subsurface)  by  migration.  Now,  we  will  show  how  the  geometry  of  reflector 
is  determined  by  the  geometry  of  the  traveltime-offset  curve. 

We  shall  denote  by  AT  a  2-D  vector,  X  =  ( x ,  z)  .  Let  xs  be  source  position  and  xT 
be  receiver  position  located  on  the  horizontal  datum  surface  L  with  y  the  midpoint 
and  h  the  half-offset: 

xs  =  y-h,  xr  =  y  +  h. 

For  any  point  X  below  the  surface,  ts(xs,  X)  and  Tr(X,xr),  respectively,  denote  trav- 
eltimes  from  xs  to  X,  and  from  X  to  xT. 

Suppose  we  know  the  total  reflection  travetime  T(y ,  h)  and  a  background  velocity 
c(x,  z).  Then,  for  each  half-offset  h ,  the  reflector  is  determined  such  that 


Tt(xs,  X)  +  Tr(X,  xr)  =  T(y,  h), 


(1) 


drs  drr  _  dT 
dy  +  dy  dy' 


(2) 


where  X  =  (x,  z)  is  the  point  on  the  reflector.  Equations  (1)  and  (2)  show:  if  we  know 
the  traveltime-offset  curve  T(y,h)  (therefore,  dT/dy)  for  some  h,  we  can  compute  a 
imaged  depth  z  for  a  fixed  location  x.  If  c(x,  z)  equals  the  true  velocity,  then  the 
imaged  depth  z  is  independent  of  offset  h;  otherwise,  for  wrong  background  velocity, 
z  varies  with  offset  h.  Consequently,  the  imaged  depth  provides  us  information  on 
velocity. 


Equations  (1)  and  (2)  are  for  the  common  offset  case, 
reflector  equation 

t,(x„  X)  +  Tr(X,  xr)  =  T(y,h), 
drr  =  OT 
dxr  dxr  ’ 


Similarly,  we  have  the 

(3) 

(4) 


for  the  common  shot  case;  the  equation 


T,(x„X)  +  Tr(X,xr)  =  T(y,h), 


(5) 


dr,  dT 
dx,  dx9  ’ 

for  the  common  receiver  case;  the  reflector  equation 


(6) 


T»(x»,X)  +  Tr(X,xr)  =  T(y,h),  (7) 

drr  dr,  _  dT 

dh  dh  dh'  K  ’ 

for  the  common  midpoint  case.  Among  these  cases,  only  common  offset  has  a  sym¬ 
metric  imaged-depth  function,  which  allows  for  the  possibility  of  approximating  the 
imaged-depth  function  by  a  hyperbola. 
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Specifically,  when  the  background  velocity  is  the  constant,  c,  then 

Ta(X,  Xs)  —  pt/c,  Tr(X,Xr)=Pr/c, 

where  _  _ 

pa  =  \j{xa  -  x)2  +  22,  Pr  =  \J{xr  -  x)2  +  Z2. 

For  this  case,  equations  (1)  and  (2)  are  simplied  to 

p»  +  Pr  =  cT{y,h),  (9) 

do) 

dy  dy  dy 

Next,  we  will  study  the  quantitative  properies  of  z(h)  when  c(x,  z)  differs  from  the 
true  velocity. 


RESIDUAL  MOVEOUT 

When  the  background  velocity  c(x,z)  differs  from  the  true  velocity,  there  is  a 
deviation  between  the  imaged  depths  of  different  offsets;  i.e.,  a  residual  moveout  from 
a  horizontal  alignment  is  observed.  We  expect  that  this  residual  moveout  can  be  used 
to  measure  the  error  of  the  velocity  and  be  independent  of  the  dip  of  a  reflector. 

Constant  Velocity  and  Horizontal  Reflector 

Suppose  that  the  true  velocity  and  background  velocity  are  constants  and  the 
reflector  is  horizontal.  In  this  special  situation, 

T2(y,  h)  =  T2(y,  0)  +  4/i2/ v2 , 

P»-  Pr  =  y/(z2  +  h2)  =  |r(y,  h ), 

so  that 

z\h)  =  jTHy,h)  -h?  —  jT°(y,  0)  +  (c2/^  -  1  )h3  =  z2(0)  +  (<?/v2  -  1  )h2. 

That  is,  the  moveout  is  the  hyperbola 

z2(h)  =  z2(0)  +  (c2/v2-l)h2.  (11) 

Equation  (11)  shows  that  the  residual  moveout  is  the  exact  hyperbola  for  constant 
velocity  v  and  horizontal  reflector.  This  result  is  parallel  to  that  of  moveout  in 
unmigrated  data. 

General  Case 

For  general  velocity  or  arbitrary  reflector,  one  should  not  expect  a  simple  expres¬ 
sion  such  as  (11).  Instead  of  that,  we  consider  the  asymptotic  expression  under  the 
assumption  of  small  offset. 
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First,  because  z  is  a  symmetric  function  of  h,  we  get 


dz  . 

dii  u-0=  °' 

This  implies  the  Taylor  series  expansion 


d3z  . 

dh 3  lta*=  0l  elc- 


1  tfiz2 

z*(h)  =  z2(0)  +  i-jjj-  |tao  h2  +  0(h*).  (12) 

Now  let  us  try  to  estimate  the  second  derivative, 

dV  . 

dh 2  1/1=0  ' 

For  fixed  x,  the  midpoint  y  and  the  imaged  depth  z  are  functions  of  offset  h. 
differentiating  equation  (1)  with  respect  to  h,  we  have 

\drs  drA  dy  \dra  drrl  \dj±  diA  dz_  _  dTdy  dT 

I  “T  I*  T  T  ft,  +  U  +  I  J,  11  »  •  V^/ 


0y  dyl  dh  I  dh  dh  I  [dz  dzj  dh  dydh  dh 


Using  (2),  we  get 


dr,  drr  dz  _  dT  dra  drT 
dz  dz  dh  dh  dh  dh 


Notice  that  y  is  symmetric  in  h,  so 


ii— • 


Lemma.  Suppose  that  /(x,  y,  z)  is  a  smooth  enough  function.  For  any  function 
y  =  y(x)  and  z  =  z(x),  if 

df  l*=0=  °»  Tx  |j:=0=  °’ 


^/(x,y(x),z(x))  U=o=  ^  U=o  • 


Proof.  By  the  definition  of  the  total  differential, 


dfdy 


dydx 


±  U=  £  U  +  |£5I  U-o  +  l^sl  I—  £  U 


d/dz 


This  completes  the  proof. 
Using  this  lemma,  and 


ss  |**°“  °- 


dz  . 

dh  |fc=0=  °* 
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we  get 


/ar 

dr3  drr 

\  |  _  a2r 

a2r,  52rr 

[dh 

dh  +  dh 

J  U=0-  a(l2  U=0 

dh2  +  dh2 

Also,  from  |/,=o=  0,  we  have 


d_ 

dh 


'(dr,  drr  dz 
dz  dz  dh 


9rt  drr  <Pz 
[dz  dz’dh 2 


(15) 


(16) 


Differentiating  equation  (14)  with  respect  to  h,  and  using  (15)  and  (16),  we  set  up 
the  following  equation  for  dPz/dh2  |/,=0 


(dr,  drr  <Pz 
K  dz  dz  ’  dh 2 


i  _ d2T 

k=°-  dh2 


d2T,  d ^7V 

dh2  +  dh2 


|/»=o  • 


(17) 


Equation  (17)  holds  for  any  velocity  function  v(x,  z),  any  background  velocity  c(x,  z), 
and  an  arbitrary  reflector.  Now  we  will  simplify  equation  (17)  for  the  laterally  invari¬ 
ant  background  velocity. 

Suppose  the  background  velocity  is  the  constant,  c.  Then 

ts  =  yjz2  +  (x-  y  +  h)2/c ,  rr  =  \jz2  +  (x  -  y  -  h)2/c. 


After  calculation, 


where 
We  have 


dr ,  drr  z 

r\  I  h=Q  ”  r\  /l=0“  * 

dz  dz  cp 

d2rs  _  dPTr  _  z2 
dh 2  'h=0~  dh2  'h=0~  cp3!2' 


p  =  yj  z2  +  (x  —  y)2. 


dr3  drr 
dz  dz 


d2rt  d2rr  2_  _  2(x  -  y )2 

dh2  dh2  /l=0  cp3  cp  cp3 


Using  these  results,  equation  (17)  is  simplified  by 


Furthermore, 


&£*M  _  i  ,  2(*  -  y )2  2_ 

cp  dh2  ^=0  dh2'h=0+  cp3  cp ' 


2z 


d?z 

dh2 


I  =  —  1 

|  h=0  l^=0’ 


(18) 
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P  =  2T  U=o, 

( V  ~  x)  _  cdT 
p  2  dy  U=0  ' 

Thus,  we  obtain  the  result 


or 


where  r  is  the  migration  time. 


(19) 


(20) 


From  equation  (19),  we  conclude  that  after  migration,  the  the  main  part  of  the 
residual  moveout  is  determined  by  the  background  velocity  c,  traveltime  T,  curvature 
of  unmigrated  data  d2T/dh2,  and  slope  of  unmigrated  data  dT/dy.  This  is  true  for  a 
constant  background  velocity  c,  any  velocity  function  v(x,  z),  and  arbitrary  reflector. 
In  addition,  the  stacking  velocity  is  defined  by 


1  _  1 
[ca<*]2  4 


(21) 


which  can  be  directly  estimated  from  the  residual  moveout.  Compared  to  NMO,  we 
have  the  new  term,  dT/dy.  In  the  follow  examples,  one  will  see  this  term  removes 
the  dip  effect  in  RMO  velocity  analysis. 

For  a  small  gradient  velocity  c(z)  or  a  small  dip  of  the  reflector,  equation  (20)  can 
be  modified  to 


cPr* 
~dh T 


where  c„„,  is  the  root-means-squared  velocity  of  c(z). 


(22) 


Examples  For  Constant  Background  Velocity 

We  use  constant  background  velocity  for  these  examples.  For  arbitrary  c(x,  z),  the 
residual  moveout  formula  will  become  too  complex  to  do  velocity  analysis.  Further¬ 
more,  prestack  migration  is  fastest  for  constant  velocity.  In  addition  to  the  normal 
moveout  formula  (11),  we  now  compute  residual  moveout  by  the  formula  (19)  and  (20) 
for  several  special  cases.  These  results  are  similar  to  normal  moveout,  except  that 
the  former  is  insensitive  to  the  dip  of  the  reflector. 
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1.  Constant  velocity  and  dipping  reflector 

Suppose  that  the  true  velocity  is  a  constant,  v ,  and  the  reflector  dip  is  9.  In  this 
case, 

T2(y,  h)  =  4 y2  sin2  9f  v 2  +  4/i2  cos2  9/v2. 

We  have, 

U=o=  4  cos2  9/v2, 
y-  U=0=  2  sin  9/v. 

Therefore, 

, a, 2 

-ftp  U= 0=  2(c2/v2  -  1).  (23) 

z2(h)  =  22(0)  +  (c2/v2  -  1  )h2  +  0(h4).  (24) 

or 

r2(/i)  =  r2(0)  +  (1/v2  -  1/02)4 h2  +  0{h4).  (25) 

Equation  (25)  shows  that  the  residual  moveout  is  independent  of  the  reflector  dip 
and  the  stacking  velocity  is  v,  when  the  true  velocity  v  is  a  constant. 


tZL 

dh* 


2.  Constant  velocity  and  diffraction  from  a  scattering  point 

Suppose  that  the  true  velocity  is  a  constant,  v,  and  a  scatterer  is  located  at  the 
point  ( xm,z* ).  In  this  case, 


T(y,h) 


We  have, 


Therefore, 


or 


\J(y  —  h  —  x‘)2  +  (zm)2  +  yj(y  +  h  -  xm)2  +  (2*)2]  / v 


Lfll  ,  _  4  (**)a 

[  a*2J  ,'*=0  u2(y-x*)2  +  (z')2’ 

/ dT\ 2  4  (y  -  x*)2 

\dy)  U=0  v2  (y  -  x*)2  +  (x*)2’ 


U=o=  2(c2/n2  -  1). 

z2(h)  =  22(0)  +  (c2/v2  -  1  )h2  +  0(h4). 
r2(h)  =  r2(0)  +  (1/v2  -  l/c2)4/i2  +  0(h4). 


(26) 

(27) 

(28) 


Again,  equation  (28)  shows  that  the  residual  moveout  is  independent  of  the  lateral 
offset  from  the  point  scatterer  and  the  stacking  velocity  is  v,  when  the  true  velocity 
v  is  a  constant. 
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3.  Laterally  invariant  velocity  and  horizontal  reflector 

Suppose  that  the  true  velocity  is  a  laterally  invariant  function,  v(z),  and  the 
reflector  is  horizontal,  with  depth  zm.  In  this  case, 

T2(y,h)  =  T2(y,  0)  +  4/i2/[tvmj(**)]2  +  0(/i4). 


We  have, 


U=0=  4/[Vrm,(0]2, 

dT  _ 

¥  U—  0. 


Therefore, 

£  '-= 2  {idpw  -  0  •  (29) 

*2(A)=zWidhr0A2+o(^  <30) 

Notice  that  z(0)  ^  z*.  That  is,  in  depth  migration,  the  imaged  depth  is  inconsistent 
with  the  desired  point  at  which  the  root-mean-square  (RMS)  velocity  is  determined 
from  the  residual  moveout.  However,  if  we  let 


r*  =  2  / 
Jo 


■**  ds_ 

•  v(sy 


tnen 

T2(',)=r2<0)+(id^-?)4'‘,+O(n  <31) 

and  r(0)  =  r*.  Therefore,  time  migration  can  give  us  the  correct  location  at  which 
the  RMS  velocity  is  determined  by  the  residual  moveout. 

Equation  (31)  shows  that  when  tl:  ?  true  velocity  is  laterally  invariant,  the  stacking 
velocity  determined  from  RMO  is  consistent  with  RMS  velocity  at  the  imaged  time. 


4.  Laterally  invariant  velocity  and  dipping  reflector 

Suppose  that  the  true  velocity  is  a  laterally  invariant  function,  v(r),  and  the 
reflector  is  a  dip,  with  angle  9.  In  this  case, 

=  l  (l-pV(<r))  1/2  da, 
d2T  4 

dh2  l**°-  jr-  V2(l  _p2v2)-3/2d(T’ 

dT  _9n 
dy  l*-°” 
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where  r*  is  the  vertical  time  at  the  reflection  point  and 

sin  6 


P  = 


v(r*)' 


Therefore, 


(a 

T\ 2 

T  dh2  +  (c 

}y) 

I  .  ,  /,r'(l-pV)-]/2fe 

h  °  fo*  v2(l  —  p2v2)~3/2dcr 


The  stacking  velocity  is  estimated  by 


r  ,2  /o’^)(l-AV))  3/2<*g 


(32) 


(33) 


Notice  that  cstk  does  not  equal  the  RMS  velocity  v 2  that  is  defined  by 

vl  =  ~:J0  v2(a)da.  (34) 

In  fact, 

[cstk]2  =  v\  +  |(w$  -  v$)p2  +  0(p4),  (35) 

where 

v\  =  —  /  v4(s)ds. 

Tm  JO 

v4  is  always  greater  than  v2  and  they  are  equal  only  if  v(z)  is  a  constant.  Therefore, 
c,tk  is  always  greater  than  v2  and  they  are  close  for  small  p  or  a  small  gradient  of 
v(z).  Furthermore,  the  imaged  time,  r(0),  is  different  from  r*,  and 

[i’2(t(0))]2  =  v\{Tm)  +  ^(v|(r‘)  -  c2)(v2{tx)  -  vf(r*))p2  +  0(p4).  (36) 


Usually,  v  is  bigger  than  v2,  so  v2(r(0))  is  bigger  than  v2(rm)  when  c  is  smaller  than 
v2(t").  If  c  is  chosen  suitably  small,  cstk  may  be  a  good  approximation  to  v2(r(0)). 

we  conclude  that  for  a  suitable  background  velocity,  the  stacking  velocity  from 
RMO  velocity  analysis  may  approximate  the  RMS  velocity  at  the  imaged  time. 


Higher  Residual  Terms 


The  residual  moveout  is  hyperbola-like  only  for  small  offset.  In  fact,  offsets  should 
not  be  too  small  so  that  we  can  have  a  high  resolution  in  velocity  analysis.  (See  Liu 
and  Bleistein.)  Therefore,  we  require  an  error  estimate  for  formulas  in  nonzero  offset. 
This  work  is  partially  implemented  in  Mathematica. 


Using  the  rule  of  differential  for  a  compound  function  in  equation  (14)  and  setting 
c  constant,  the  fourth  order  derivative  of  z  with  respect  to  h  satisfies 


dAz 


d(p»  +  Pr) 


dz 


d*T  (Py&T 
dh*  +  3  dh2  dh2dy 


U-o=  c 

(Pz  &{ps  +  Pr) 


u- 


&{p>  +  Pr) 


-6 


dh2  dh2dz 


U=o  -3 


dh* 

<Py  &{p»  +  Pr)] 


U=o 


dh 2  dh2dy 


|fc=0 


(37) 
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For  a  constant  velocity  v  and  a  dipping  reflector  with  angle  9 ,  we  can  obtain 

dV  _  24(1 /u2  -  l/c2)sin2  29 

dh4  'h~°  (2*)2  ’  (38' 

where  r  is  the  migration  time  and  zm  is  the  reflection  depth.  Therefore,  from  (38) 
and  (23),  a  more  accurate  expression  for  the  residual  moveout  is 

r'W  =  *0)  +  4  (1  -  i)  M  -  +  0,h«) 

=  r2(0)  +  4  (i  -  i)  A2  (l  -  j(p)2sin22^  +  O(V).  (39) 

For  a  constant  velocity  u  and  a  scattering  point  at  (x*,  zm), 

|  24(1 /u2  —  l/c2)sin22ff 

d/l4  ’/‘=0  (2*)2  ’ 


where 


0  =  arc  tan 


y-  x 


By  the  way,  for  the  unmigrated  data,  the  higher  residual  term  is 

d*T*  ,  24  sin2  29  cos2  9  ,  j x 

1^4"  1*=°“  V2(2.)2  *  (41) 

Equations  (38)  and  (40)  show  that  when  the  true  velocity  is  a  constant,  small  higher- 
residual  terms  of  moveout  are  obtained  for  a  closed  background  velocity,  small  ratio 
of  offset  to  the  imaged  depth,  and  the  dipping  angle  that  is  near  0  or  90  degree. 

For  laterally  invariant  velocity  and  a  horizontal  reflector  with  depth  zm, 


d*T2  ,  24(^2  -  vi) 

dh4 


where  to  is  the  zero-offset  time,  and 


(42) 


V2  =  r  f  v2(s)ds, 

1 0  Jo 

v\  =  —  f  vA(s)ds. 
to  Jo 

The  result  in  (42)  is  the  same  as  in  the  unmigrated  data. 

Equation  (42)  shows  that  when  the  true  velocity  is  laterally  invariant,  a  small 
higher-residual  term  is  obtained  for  small  gradient  of  the  velocity. 
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VELOCITY  CONVERSION 

In  the  previous  chapter,  we  show  how  stacking  velocity  is  obtained  by  moveout 
methods.  However,  an  interval  velocity  function  is  required  for  migration  and  inver¬ 
sion.  Usually,  we  assume  that  the  stacking  velocity  equals  the  RMS  velocity.  Then 
the  interval  velocity  is  computed  by  Dix  equation  or  other  algorithms.  Unfortunately, 
this  assumption  fails  when  the  velocity  has  a  lateral  anomaly.  Now  we  will  give  an 
equation  to  solve  for  the  laterally  varying  interval  velocity.  This  equation  requires 
the  small  lateral  variation  and  horizontal  reflector. 

Suppose  the  true  slowness  w(x,  z)  can  be  written  as 

w(x,z)  =  w(z)(l  +  a(x,z)),  (43) 

where  w(z )  is  a  reference  slowness  and  a(x,  z)  is  a  small  perturbation.  We  obtain  the 
equation 


6ws(y,z)  = 


T0v3{z) 


+  a  (1  + 


£ V). 


T  <44> 


where  6w,  is  the  anomaly  of  the  stacking  slowness,  To  is  the  zero  offset  time,  and  v3  is 
the  RMS  velocity  from  the  referenced  velocity.  When  a  and  w  are  depth  independent, 
the  result  in  equation  (44)  is  the  same  as  that  of  Lynn  and  Claerbout.  Equation  (44) 
shows  that  the  second  derivative  of  a  determines  the  anomaly  of  the  stacking  velocity; 
the  anomaly  of  the  stacking  velocity  at  a  depth  results  from  the  anomaly  of  the 
interval  velocity  above  this  depth.  Furthermore,  //  vds  increases  as  a  decreases,  so 
the  anomaly  of  the  interval  velocity  near  the  surface  has  the  largest  effect. 

Applying  Fourier  transform,  with  respect  to  y,  to  equation  (44),  we  obtain 


6wa(kx,z)  = 


Tvt{z) 


£  i*..')  [-*:  (4jf ) 


♦'■•S3' 


dc 

(45) 

v 


Equation  (45)  is  a  First-kind  Volterra  integral  equation  that  is  ill-posed.  Therefore, 
the  recursive  alogrithm  for  equation  (45)  is  unstable.  To  obtain  a  stable  solution,  we 
may  apply  the  damping  least-squared  method  to  equation  (45). 


SYNTHETIC  DATA  EXAMPLE 

To  test  the  residual  moveout  method,  we  applied  this  method  to  synthetic  seis¬ 
mograms  computed  for  a  subsurface  model  in  which  velocity  increases  linearly  with 
depth  z ,  according  to  v(z)  =  1.5-f  0.8z  km/s.  The  model,  shown  in  Figure  1,  consists 
of  five  reflectors,  each  with  a  dipping  and  horizontal  segment.  Dips  for  the  dipping 
segments  range  from  30  .o  90  degrees  in  15-degree  increments.  The  seismograms 
contain  10  offsets,  ranging  from  100  m  to  1900  m  in  200  m  increments.  Because  of 
dipping  reflectors  and  depth  dependent  velocity,  the  stacking  velocity  in  equation  (33) 
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and  the  RMS  velocity  are  not  same  but  close  each  other.  The  error  between  the  both, 
shown  in  Figure  2,  increses  with  depth  and  dip. 

After  prestack  migration  with  the  constant  velocity,  c  =  1.5  km/s,  one  of  the 
CDP  gathers  is  plotted  in  Figure  3.  Because  the  background  velocity  is  lower  than 
the  true  velocity,  all  event  locations  increase  with  offset.  The  velocity  scan  for  this 
CDP  is  plotted  in  Figure  4.  Unlike  the  velocity  scans  in  NMO,  velocity  peaks  here 
are  single- valued.  After  residual  moveout,  all  events  are  corrected  to  horizontal  ones 
(in  Figure  5).  The  stacking  result  is  shown  in  Figure  6,  which  is  equivalent  to  the 
poststack  migration  with  the  constant  velocity,  c.  By  using  the  interval  velocity 
converted  from  the  stacking  velocity,  poststack  residual  migration  gives  the  correct 
reflector  positions  (in  Figure  7). 


Model 


Fig.  1.  Subsurface  model  used  to  generate  synthetic  seismic  traces. 


CONCLUSION 

Velocity  analysis  by  prestack  migration  can  handle  dipping  reflectors.  Conven¬ 
tional  approaches  use  iteration,  which  results  in  larger  computation  than  the  method 
proposed  in  this  paper.  Using  the  relationship  between  the  residual  moveout  and 
the  error  in  the  background  velocity,  we  can  estimate  directly  the  true  velocity  with¬ 
out  iteration.  Furthermore,  stacking  after  the  residual  moveout  provides  a  partial 
migration  output.  With  this  output,  a  residual  post  migration  should  yield  a  more 
accurate  structural  image.  Using  a  more  general  background  velocity,  we  may  handle 
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Dip  (degree) 

0  30  60  90 


Contours  of  Velocity  Error 


Fig.  2.  The  relative  error  between  the  stacking  velocity  and  RMS  velocity.  The 
difference  of  contours  is  0.002.  The  arrow  direction  indicates  increase  of  the  error. 


Common  location  images  (CDP  =  161 ) 


FIG.  3.  One  of  the  CDP  gathers  after  the  migration  with  the  constant  velocity. 
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Velocity  (m/sec) 

1500  2500  3500 


Velocity  Spectrum  CDP161 


Fig.  4.  Velocity  analysis  for  the  CDP  gather  in  Fig.  2. 


Offset  (km) 


0.5  1.0  1.5 


Fig.  5.  Residual  moveout  for  the  CDP  gather  in  Fig.  2. 
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Midpoint  (km) 

-10  12  3 


Stack  of  CDPs  1  to  301 


Fig.  6.  Stacking  for  the  ten  offsets. 

Midpoint(km) 


Residual  Migration 


Fig.  7.  Poststack  residual  time  migration. 
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the  lateral  variation  of  velocity,  but  it  will  require  a  complicated  algorithm.  The 
perturbation  method  here  is  suggested  to  handle  the  lateral  variation  of  velocity. 
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APPENDIX  A:  DERIVATION  OF  EQUATIONS  (35)  AND  (36) 
From  Taylor’s  expansion, 


(1  -  pV(<r ))  3/2  =  1  +  - v2{a)p2  +  0(p 4), 


so  that 


(A-l) 


JQ  v2(<r)(l  -  p2v2(a))~3/2da  =  jf  v2(a )d<r  +  ^P2JQ  v\a)da  =  rm(vj  +  ^v*p2), 

JQ  (1  -  p2v2(a))-3/2d<r  =  r*(l  +  ^v\p2). 

Using  these  results  and  the  definition  in  (33),  we  obtain 

fc  ,J2  -  r*(t;2  +  3/2  v\p2)  _  2  3,  4  4  2  .  4- 

lCj,fcJ  "  r*(  1  +  3/2  ^Ip2)  2  +  2(  4  ^  + 

This  completes  the  proof  of  equation  (35). 

From  the  Taylor’s  expansion, 

dv2 

”2 (t)  =  vJ(t')  +  (r  -  |r=r.  +0((t  -  r")2) 

=  <’!(’■•)  +  (r/f  -  1)(«V)  -  ^(t*))  +  0((r  -  r')2), 

where  we  use  the  fact  that 

dv2 

-Jl  |r„.=  („V)  _  „2(r‘))/r\ 

From 


(A-2) 


and 


we  have 


r(0)  =  r(»,0)^l  -c2?2  =  r(y,0)(l  -  ic2?2)  +  0(p4), 

^ts/.O)  =  J0  (1  -pV(<r)r'/2rf<7  =  T-(l  +  ivfp2)  +  0(p4), 

r(0)  =  p-(l  +  i(»|(r-)  -  c2)?2)  +  0(p4). 

Subsitituting  the  above  formula  into  equation  (A-2)  and  setting  r  =  r(0),  we  obtain 

uf(r(0))  =  v2(rm)  +  ^(u| (rm)  -  c2)(v2(rm)  -  v\(rm))p2  +  0(p4).  (A-3) 

This  completes  the  proof  of  equation  (36). 
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APPENDIX  B:  DERIVATION  OF  EQUATION  (44) 

Suppose  reflectors  are  horizontal  and  the  true  velocity  v(x,  z)  can  be  written  by 


v(x,z)  v(z) 


(1  +  a(x,z)), 


(A-4) 


where  v(z)  is  a  reference  velocity  and  a(x,z)  is  a  small  perturbation.  Under  the 
assumption  of  the  small  perturbation,  we  can  calculate  the  two-way  traveltime  by 


TiKy)=TM+2j;^lda 


(A-5) 


where  9  is  the  angle  of  the  ray  path  from  the  vertical,  (f ,  a)  is  a  point  on  the  ray 
path,  and 

£(<7 )  =  y  +  J  tan  Bds.  (A-6) 

Therefore, 

d2T{h,y)  ,  d2T(h,y)  ,  ,  da 

~dh?  u=0=  l*-°  +2Jo  d&  [l^r)  1,1=0  T  (A'° 

From 

d2  fa(Z,a)\  d2a  1  da  9  (  1  \  ,  ,,  .  d2  (  1  \  ..  0. 

ah2\  cos 9  )  ~  dKicos9+2dhdh\cosd)+a^'<T)dh?\cos0j'  (  "  ] 

and 

A  (-!_)  |4=0=  o, 

ah  Vcos 9) 
cos  9  1^0=  1, 


we  have 


w  i‘-=  w  I-  +°<*  '>&  (sb)  I-  •  (A-9> 


Notice  that  £  is  a  function  of  h,  so  that 


_  a2a  /df\2  da  d2^ 
~  d?  \dh)  +  d(dh2' 


From  (A-6),  £  -  y  is  an  odd  function  of  /»;  hence 


v  s  I  _ 

dh 2 1,1=0 

This  result  and  formula  (A-10)  give 


d*(t  -  y) 


Ikm  0=  0. 


d2a  .  _  d?a 

ah2  dy2 


W  U=0- 


(A-10) 


(A-ll) 
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Here  we  use  the  fact  that  £  =  y  at  h  =  0.  To  do  the  further  calculation,  we  introduce 
the  slope  parameter 

sin# 

P  = 


via)' 


which  is  independent  of  a.  Then 


£(<r)  =  y+  f  ,  VP  ds , 

Ja  \Jl  -  (vp)2 


_  r  vp 

J°  \Jl-  (vp)2 


ds. 


These  formulas  imply  that 


d£  _  rz  v 

dp  Ja  (l  -  (up)2)3/2 

dh  _  [*  v 

dp  Jo  (1  —  (up)2)3/2 


ds, 

ds. 


dZ  |  _  d£  .dh.  _  fi  vds 

dh  ^=0  dp  '  dp  I'1"0  fi  vds' 


From 


we  have 


1 


1 


d2 
dh 2 


cos*  ^l-(up)2’ 

(cos#)  ”  dh2  (^1  -  (up)2) 


__  d?_  (  1  \  /  dp  \ 2  d_  (  1  \  d2p 

~  dP*  \y/l-(vp)2)  \Sh)  dP  \  y/l -{vp)2J  dh2' 


Again,  p  is  an  odd  function  of  h,  so  that 

92P  |  _  Q 

ai? Ita0_  °- 

This  result,  formulas  (A-12),  and  (A-14)  give 

fl*  /  1  \  i  &  ( i  _i _ gy) 

|ta0  dr?{fizj?pp)  (MW 

By  formulas  (A-9),  (A-ll),  (A-13)  and  (A-15),  we  obtain 


(A-12) 

(A-13) 


(A-14) 


(A-15) 


dh2\cos9  )  |a=0  dy 2  (finds)  +  (v,  \ finds)*'  (A'16) 
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Substituting  this  formula  into  equation  (A-7),  we  have 


d2T(h,y)  d2T(h,y) 

dh 2  ]h=0  Oh 2 


l  [Z  ^2q  (fgVds\2  <xv2(° 
h~°  Jo  dy2  [fjvdsj  (So  vds 


2(o)  do 
vds)2  v 


(A-17) 


The  stacking  velocity  vs  is  defined  by 


_i _ i  \T,  ,,fr(M)l . 

v](y,z)  4  T  *  a/>2  l‘"° 


(A- 18) 


By  using  formula  (A-5)  and  (A-17),  we  have 


tw]  '-»=  [fw]  <*-  +2w + 


(A-19) 


From 


Tsm  _j_ 
m  l*_0  aw 


(A-19)  becomes 


av2(o)  do 
fj  vds)2  v 


_  4  _  __8_  r*  ado_  ^  y*  / J^Vds\2  _ai 

v2s(z)  Tv2  Jo  v  Jo  dy2  {j0zvdsJ  (f0z 


£  ,  .  1  f*  ado  Tvs  f*  \dPa  ( flvds\2  av2(o)  do 

ws(y,z)  -  TijJo  —  +  —J0  [dy2  +  (JS  vds)2 j  T 


(A-20) 


where  we  use  the  relation 


atrr.fy,^)  =  — 1 _ !_  =  M£l(_2 _ M 

v.(y,z)  v,(*)  2  Vt;2(y,z)  v2  (z)J 


From 


I' vd,  =  £f>?(r). 


equation  (A-20)  becomes 


,  ,  .  1  [*  f^a  / /*  vds\ 

w‘(y’z)-  Tv.(z)h  [aj?  U.wJ 


-l-a  1  +  TjH  — . 

v2(z)  V 


(A-21) 
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